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ABSTRACT 

Recently the ARCADE 2 experiment measured the cosmic radio background 
(CRB) and found the brightness temperature of the CRB at 1.4 GHz to be 
~480 mK. Integrating the flux density from the observed 1.4 GHz radio source 
count produces a brightness temperature of ~100 mK — less than a quarter of the 
observed CRB at 1.4 GHz. Radio quiet AGN are a large fraction of the 1.4 GHz 
/iJy sources and typically host significant star formation. Thus, it is possible that 
AGN and host star formation could be responsible for some fraction of the excess 
CRB at 1.4 GHz. Here, an X-ray background population synthesis model is used 
in conjunction with empirical radio to X-ray luminosity ratios to calculate the 
AGN contribution to the CRB at 1.4 GHz including the emission from host star 
formation. It is found that AGN and host star formation contribute <9% of the 
CRB at 1.4 GHz. When all known 1.4 GHz radio source classes are considered, 
<60% of the CRB at 1.4 GHz is accounted for; therefore, it is necessary that 
either known radio sources evolve significantly at flux densities below current 
survey sensitivity limits or a new population of low flux density radio sources 
exist. 

Subject headings: galaxies: active — galaxies: Seyfert — quasars: general — 
galaxies: star burst — radio continuum: galaxies 



1. Introduction 



Extragalactic backgrounds are important tools in understanding the density and evo- 
lution of different objects with in the universe. T he 7-ray background places constraints on 
dark matter annihilation (e.g.. ICirelli et al.ll201Cll ). the X-ray background (XRB) provides a 
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census of accretion onto supermassive black holes (e.g.. lGilli et al.ll2007l ). the cosmic infrared 
back ground (CIRB) e ncodes the history of obscured star formation and galaxy evolution 



[e.g., iLutz et al.l 1201 if ), and the cosmic r aicrowave background (CMB) provides the deep- 
est insights into the early universe (e.g., iKomatsu et al.l 1201 ll ). It is known that at least 
two source populations contribute to the cosmic radio background (CRB). At fluxes > 1 
mJy radio loud active galactic nuclei (AGN), actively accreting supermassive black holes 
defined by f^l fpt 10^ w here /s is the 5 GHz flux density and fs is the B band flux densit y 
( iKellermann et al.l 1 19891 ). are the dominate source population (e.g.. IWindhorst et al.lll993l ). 
At lower flux densities, the 1.4 GHz source count steepens as star forming galaxies (SFGs) 
and radio quiet AGN begin to dominate. Recently, the ARCADE 2 experiment measured the 
CRB at frequencies between 3 and 90 GHz and found that the spectrum of the observed CRB 
suggests that the extragalactic brightness temperature at 1.4 GHz is ~480 mK (Fixsen et 
al. 2011, their eq. 6). Even when considering the large dispersion between various observed 
source counts, extrapolating the total observed source cou nt at 1.4 GHz to lower flux densi- 
ties can only account for ~1 00 mK of the observed CRB (iGervasi et al.ll2008t ISeiffert et al. 



201 ll : IVernstrom et al.ll201l[ ). The next step must be to consider phenomenological models 
of individual source classes. Here, the contribution of AGN and their hosts to the CRB is 
investigated. 

Radio quiet AGN, although less radio luminous, are much more common than radio 
loud AGN, and thus are an imp ortant class of sources in the /xJy flux density range at 1.4 



Lutz et al. 


201ol 


Yuan et al. 


2010) 



GHz (e.g., 

to host significant star formation (e.g. 
plausible that emission from radio quiet AGN and their hosts makes a significant contrib ution 
to the CRB at 1.4 GHz. Using optical AGN luminosity functions, ISingal et al.l (12010! ) find 
that radio quiet AGN should contribute a few percent of t he CRB. While optical AGN 
luminosity functions do include emission from the host gal axv (iDraper fc Ballantvnell2011bl) . 



star formation in AGN hosts tends to be obscured (e.g., iLutz et al.ll2010l ): therefore, AGN 
host star formation contributes very little to the total optical luminosity of the AGN system. 
Star formation in AGN hosts must be carefully considered in order to understand the possible 
range of CRB fractions which can be accounted for by AGN system radio emission. 

Here, a detailed accounting of the AGN contribution to the CRB is conducted with 
the inclusion of star formation in AGN host galaxies. The framework of XRB population 
synthesis modeling is combined with empirical X-ray to radio flux ratios. This allows for 
the consideration of AGN with radio flux densities well below the current survey sensitivity 
limits. Variations in the AGN contribution to the CRB due to uncertainty in the evolution 
of AGN host star formation, the AGN luminosity function, and the AGN radio to X-ray 
luminosity ratio are also explored. A A-dominated cosmology is assumed with /i = 0.7, and 
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index a, such that S oc u". 



1 — Qrn = 0.7 (ISpergel et al.ll2003[ ). Radio sources have a flux density with spectral 



Calculations 



The AGN number counts at 1.4 GHz are calculated between S„ 



10 mJy and S'„ 



= 100 nJy, in the same manner as by iBallantynd ( 120091 ). An XRB synthesis model is used 
to characterize the space density of AG N, the AGN type 2/type 1 ratio, f^^, and the fractio n 



of Compton thick (CT) AGN, fcr (e.g.. iBallantvne et al.ll2006l : lOraper fc Ballantvnell2009h 



AGN with obscuring c olumn densities Nh > lO^'^ cm~^. Either the lUeda et al.l (120031 ) or 
La Franca et al.l (120051 ) hard X-ray luminosity function (HXLF) is used to model the space 
density and evolution of AGN. The evolution of /2 is depende nt on both Ly, th e AGN 2-10 
keV luminosity, and redshift, z, in the manner determined by iBallantyne et al.l (120061 ). The 
fcT is det ermined by fitting the peak of the XRB at ^ 30 keV and the local space density of 
CT AGN dPraper fc Ballantvnel]201ol : iBallantvne et al.lboilh . Both the scenario where CT 



AGN evolve like less highly obscured type 2 AGN (10 < A'^^^ < 10 c m ) and the scenario 



wher e CT AGN are in specific evolutionary states are considered (see iDraper fc Ballantyne 
2010h . 



Following iTerashima fc WilsonI (j2003[ ). the ratio of 5 GHz radio luminosity to Lx for 
radio quiet AGN, defined as Rx = \og{uLy{hGHz) / Lx) is characterized as 



R 



X 



-0. 67 log + 23.67 

-5 

log Lx - 49 
-4 



41.5 < logLx < 43 

43 < log Lx < 44 

44 < log Lx < 45 
logLx > 45. 



of Panessa et al 



-0.70 

Ballantynd (l2009l ) showed that the radio quiet AGN Rx relation 



The AGN 5 GHz lu mi nosity is co n verte d to 1.4 GHz using a spectral index of a 
jKukula et al.lll998h. 



(2007) is also consistent with the current observed AGN radio number 



counts. The lPanessa et al.l (120071 ) radio quiet Rx, based on observations of local Seyferts, is 
nearly constant with luminosity. Radio loud AGN X-ray luminosities are converted to radio 
luminosities using the relation 

Rx ~ -2.0 (2) 



^If the a distribution for radio quiet AGN observed by Kukula et al. ( IQQSl ) is approximated as Gaussian 
with a — -0.7 and a — 0.6, the radio quiet AGN contribution to the brightness temperature increases by < 
1 mK. 
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( Terashima &: Wilson 2003 ). For the radio loud AGN a spectral index of a = -0.^ is used 
to convert from 5 GHz to 1.4 GHz luminosity. 

The HXLF contains both radio quiet and radio loud AGN. In order to separa te the radio 
quiet AGN and the radio loud AGN, a radio loud fraction, f^n, must be defined. iBallantyne 
(120091 ) found that radio number counts between 1 and 10 mJy were well fit assuming an 
exponential increase from /rl = 0.0175 at logLx = 41.5 to f^L = 0.10 at log Ly > 46.0, at 



all z. This is consistent with the overall radio loud fraction being 0.1-0.2 (e.g.. iDonoso et al. 
2009h . 



It is well documented that ga laxies hosting AGN also tend to host on-going star forma- 
tion (jPraper fc Ballantyndl2011al . and references therein), thus the radio flux density due to 
star formation within AGN host gal axies must be considered for an accurate calculation of 
the AGN radio number counts. The iBelll (120031 ) calibration is used to convert a star forma- 
tion rate (SFR) into a 1.4 GHz radio luminosity assuming that the infrared-radio correlation 
(IRC) do es not evolve wi th z. For the radio spectrum of star forming regions, a = -0.8 is 
assumed (lYun et al.ll200ll ). T he 1.4 GHz r adio d ifferential number counts, dN/dS, are then 
calculated as in Equation 2 of iBallantynd (120091 ). 



Once dN/dS is known, calculation of the brightness temperature, T(z/), is straight 
forward. The intensity /(i/) is related to dN/dS such that 



The Rayleigh- Jeans approximation gives 



dN 

Is 



(z/) ■ SdS. 



(3) 



T(.)^/(.)|, 



(4) 



where A = 21 cm is the wavelength and k is the Boltzmann constant. Thus, by applying 
the machinery of the XRB and using empirically determined conversions between X-ray and 
radio luminosities, the AGN contribution to the 1.4 GHz CRB is calculated. 

With the basic model ir i place, the c ontrib utio n of AGN system s to th e CRB is investi- 
gated. The results from th e lUeda et al.l (12003 ) and iLa Franca et al.l (120051 ) HXLFs and the 
Terashima fc Wilson! ( 120031 ) and lPanessa et al.l (120071 ) radio quiet Rx relations are compared. 



Constant and evolving star formation laws are explored. Both an evolving and non-evolving 



^If the a distribution for radio loud AGN observed bv iLin fc Mohn (|2007l ) is approximated as Gaussian 
with a = -0.8 and a — 0.6, the radio loud AGN contribution to the brightness temperature decreases by < 
0.5 mK. 
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model of fcT are considered. Analyzing such a large portion of the parameter space allows 
the maximum possible contribution to the CRB from AGN to be determined. The resulting 
models are then constrained by the observed 1.4 GHz number counts. 



Results 



3.1. Bare AGN 



First, the host star formation is ignored and the contribution to the CRB from bare 
AGNs is considered and shown in the upper p a rt of Tabl e [H The listed y^^^ v alues refer 
to the 15 AGN data points from ISeymour et al.l (|2008f) andlPadovani et al.l ( 120111 ) shown as 
the red, blue, and cyan data in Figure [TJ The iPanessa et al.l (120071) Rx p redicts an AGN 
bri ghtness temperature ~5 mK higher than the iTerashima fc Wilson ( 2003 ) Rx- Similarly, 
the iLa Franca et al.l (120051 ) HXLF produces an A GN brightness temp erature ~5 mK higher 



than the brightness temperature predicted by the lUeda et al.l (120031 ) HXLF. The blue lines 
in Figure [1] show the Euclidean normalized 1.4 GHz differential number cou nts for the bare 
AGN model w i th the lowest xled^ which assumes the iLa Franca et al.l (120051 ) HXLF and the 
Panessa et al.l (120071 ) Rx- Thus, bare AG N appear to only c ontribute ~4% of the CRB at 
1.4 GHz, in agreement with the results of ISingal et al.l ( 120101 ). 



The minimum X-ray luminosity of the XRB framework is L™" = 10^^'^ erg s~^, which 
corresponds to a minimum radio loud AGN 1.4 GHz radio luminosity of 1.8 x 10^^ W 
Hz~^. In the VLA-COSMOS survey, ra dio loud AGN were observed with 1.4 GHz radio 
luminosities down to ~5 x 10^^ W Hz~^ (ISmolcic et al.ll2009l) . Therefore, the 1.4 GHz radio 
loud AGN luminosity function of ISmolcic et al. is used to assess the importance of 

the low luminosity radio loud AGN. The radio loud AGN contribution to the CRB using 
both the pure densit y evolution (PDE) and pure luminosity evolution (PLE) versions of the 
Smolcic et al.l ( 120091 ) luminosity function are shown in the fifth and sixth rows of Table [TJ 
The dot-dashed line in Figure P show' s the radio loud contribution to the differential number 
counts using the ISmolcic et al.l (120091 ) PDE lumin osity function. As radio lou d AGN are X- 
ray luminous for only a fraction of their lifet ime (iDraper &: Ballantvnell2009l ). it is possible 
to assume that all X-ray selected AGN in the lUeda et al.l (120031 ) HXLF are radio quiet AGN. 
In this case, bare AGN contribute ~6-8% of the CRB at 1.4 GHz, depending on whether 
the PLE or PDE evolution is used for the radio loud luminosity function. 
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3.2. Accounting for Star Formation 



It is known that AGN hosts have star formation (jPraper fc Ballantynell2011a( ): thus, we 
now consi der the effe c t of A GN host star formation on the AGN and host contribution to 
the CRB. iBallantvnel fcoooh found that if AGN have either SFR « 2-3 Mq yr^^ or 

SFR ^ 0.25(1 + zy-^^hgLx - 40)^-^ Mq yr-\ (5) 

the predicted number counts are in good agreement with observations, as shown by the black 
lines in Figure [TJ However, as shown in the bottom half of Table [T], adding star formation to 
AGN hosts only increases the contribution of AGN and their hosts to the CRB by <10 mK. 



We also consider the evolving CT AGN model of iDraper fc Ballantynd (120101 ). in which 
CT AGN have a physically motivated Eddington ratio distribution. This model is in agree- 
ment with the XRB, the local s pace density of CT AGN, an d 2; ~ 2 mid-infrared estimations 
of the CT AGN space density (IDraper &: Ballantynd I2OIOI ). The third to last row in Table 
[1] shows the AGN contrib ution to the CRB using the evolving far model and the evolving 
SFR of iBallantynd (|2009[ ). The evolving model is in better agreement with the observed 
number counts and results in a higher brightness temperature than the non-evolving model. 

In order to evaluate the maximum possible contribution of AGN and their hosts to the 
CRB, it is assumed that AGN systems dominate the 1.4 GHz radio number counts down 
to the 10 jUJy flux dens i ty lev el. For this scenario, the radio loud AGN are accounted for 
usi ng the ISmolcic et al.l ( l2009l ) PDF luminosity function, and it is assumed that all AGN in 
the lUeda et al.l (120031 ) HXLF are radio quiet. The evolving fcT model is used and the radio 
quiet AGN hosts are assigned the SFR described in Equation [51 This maximal scenario is 
shown in Figure |2] and summarized in the last row of Table [TJ Thus, it is clear that AGN 
contribute to the CRB; however, the maximum contribution to the CRB from AGN and 
their hosts, which is in reasonable agreement with the observed source count, as shown in 
Figure [H is ~42 mK, or ~9%, of the CRB at 1.4 GHz. However, in the best fit model, as 
assessed by the minimum xledi AGN and their hosts contribute only ~28 mK, or ~6%, of 
the CRB at 1.4 GHz. 



4. Discussion and Summary 

We have shown that AGN systems contribute <9% of the CRB at 1.4 GHz. Thus the 
contribution of other source classes to the CRB must be considered. In Table [2] the known 
dominant source classes of the CRB are listed with the expected contribution of each source 
class. 
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It is known that star forming galaxies (SFGs) are an important source class in the AtJy 



and n Jy flux density ranges. Integrating the number counts predicted by the iPadovani et al. 



(I2OIII ) 1.4 GHz SFG luminosity function suggests that SFG only contribute ~3% of the CRB 
at 1.4 GHz. However, current surveys do not reach radio flux densities < 10 yuJy; thus the 
best estimates of the contribution of SFGs to the CRB, which do not require extrapolation 
of the observed source counts, are computed using the IRC and the CIRB. If the IRC does 
not evolve with z and maintains the loc ally observed correla tion at high 2;, then SFGs only 
contribute ~9% of the CRB at 1.4 GHz (IPonente et al.ll2011l) . It is unclear whether the IRC 



evolves with z or not, however, Ponente et al 



J201lh find that if the IRC evolves with z m 
a manner which is consistent with current observations, SFGs still only contribute ~14% 
of the CRB at 1.4 GHz. It is possible that low power, high redshift AGN systems, which 
observationally wo uld be difficult to distinguish from SFGs, may increase this percentage 
Jsingal et allboioh . 



As AGN, based on the calculation presented here, and SFGs, based on estimations from 
the IRC, can account for only ~25% of the CRB at 1.4 GHz, additional low flux density source 
classes must be considered. One possible source class is radio supernovae (RSNe), however, 
these objects are short-lived an d not com r non e nough to make a significant contribution 
to the CRB (ISingal et al.l 120101 ). iPadovanil (120 111 ) suggests that low radio power ellipticals 
(LRPEs) and dwarf galaxies may be important source classes at flux densities currently below 
survey sensitivities. LRPEs are predicted to have similar surface densities as AGN, and thus 



are not a dominate contributor to the CRB. Dwar 



gala xies are expected to dominate the 



source count at nJy flux densities, but IPadovanil (120111) pred i cts th at dwarf galaxies will 
only contribute <8% of the CRB. Furthermore, ISeiffert et al.l (1201 ll ) point out that for a 



population of sub-mJy point sources to account for the difference between the observed 
CRB and the brightness temperature predicted from the observed source count, the surface 
density of t his population must e xceed the surface density of galaxies in the Hubble Ultra 
Deep Field (IBeckwith et al.ll2006l ). Therefore, it is likely that the unaccounted for CRB is 
due to the evolution of known sources or an evolution of processes, such as those responsible 
for the IRC, and not due to an unknown source class. 

Despite the current lack of a dominant c andidate poin t sour ce class, it is unlikely the 1.4 
GHz CRB is dominated by diffuse sources. ISingal et al.l (I2OIOI ) considers the contribution 
from low-surface-brightness sources. Sources which may emit a large flux over an extended 
area may ha v e low enough surface brightness to be missed by current survey source counts. 
Singal et al.l (120101 ) flnds that the low-surface-brightness sources can at most contribute 
~25% of the CRB. Furthermore, the constraints placed on the CRB from the XRB and 
7-ray background are inconsistent with the CRB being dominated by diffuse emission from 
a population of relativistic electrons within the intracluster medium or intergalactic medium 
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(Sineal et al. 


2010; 


Lack] 


2011) 



to fully account for the CRB at 1.4 GHz. 

In conclusion, a detailed accounting of the AGN and host star formation contribution to 
the CRB has been carried out using the constraints of XRB population synthesis modeling 
and empirical X-ray to radio AGN luminosity ratios. This study has shown that AGN and 
host star formation contribute <9% of the CRB at 1.4 GHz. In a best case scenario, all 
sources currently considered in the literature contribute <60% of the CRB at 1.4 GHz. 
The unaccounted for radio emission must come from either a new low flux density radio 
population or from a known radio population with a strong evolution of the radio emission 
for low flux density objects. Therefore, it is necessary that radio surveys continue to explore 
the radio source count down to sub-/iJy levels in order to resolve the source class responsible 
for the unaccounted for CRB flux density at 1.4 GHz. 

This work was supported by NSF award AST 1008067. The authors thank the anony- 
mous referee for helpful comments. 
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Table 1. Contribution of AGN to the 1.4 GHz CRB 



HXLF Star Formation Law Ttotai (K) T^^ (K) xl 



Bare AGN 



U03 


0.013 


0.011 


8.3 


U03^ 


0.019 


0.011 


1.0 


LF05 


0.018 


0.015 


4.1 


LFOS'^ 


0.023 


0.015 


0.6 


S09 PDF + U03^ 


0.033 


0.031 


7.6 


S09 PLF + U03^ 


0.023 


0.021 


1.7 



AGN + Star Formation 



U03 SFR = 3 Mq yr-i 

LF05'' SFR = 2 Mq yr^^ 
U03 Eq. E 

LF05 Eq. 

U03^ Eq. E 

mS" Eq. [5p 

S09 PDF + U03^''= Eq. E 



0.021 


0.011 


1.5 


0.028 


0.015 


1.3 


0.019 


0.011 


1.9 


0.028 


0.015 


0.4 


0.024 


0.014 


0.9 


0.029 


0.014 


5.5 


0.042 


0.031 


19 



^Radio quiet AGN Lx is converted to radio luminosity using the iPanessa et al.l (|2007| ) 
conversion. For all other models the radio quiet Lx is converted to radio luminosity using 
Equation [TJ 



Radio loud AGN are accounted for using the lSmolcic et al.l (120091 ) 1.4 GHz radio luminos- 
ity function for radio loud AGN which is evolved using either pure density evolu tion (PDF) 
or pur e luminosity evolution (PLF). The radio quiet AGN are modeled using the lUeda et al. 
(120031 ) HXLF assuming that all X-ray AGN are radio quiet. 



^Uses the evolving fcr niodel of iDraper fc Ballantynd (120101 ) 
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'^Star formation follows the functional form of equation [5] with a different normalization 
factor. Compton thin AGN have normalization factor 0.3 while low Eddington ratio CT 
AGN have normalization factor 0.05 and high Eddington ratio CT AGN have normalization 
factor 4.0 



Note. — U03 refers tolUeda et all (120031 ). LF05 refers to lLa Franca et all (120051 ). and S09 
refers to Smolcic et al. ( 20091 ). Tr^r is the brightness temperature of radio loud AGN. 
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Table 2. Contributions of various sources to the 1.4 GHz CRB 



Brightness Temperature (K) Reference 



Total Measured CRB 



0.48 ± 0.07 



Fixsen et al. (2mi) 



AGN 
AGN+SF 

(1) Max AGN+SF 

SFG 

SFG non-evolving IRC 

(2) SFG evolving IRC 

(3) RSNe 

(4) LRPEs 

(5) Dwarf SFGs 

(6) Low-surface-brightness Sources 
Total (1-6) 



0.018 
0.025 
0.042 

0.015^ 
0.040^ 
0.063" 

< 0.00017 
O.OIO'^ 
0.038 
0.120 

0.27 



this work 
this work 
this work 

this work 

Ponente et al. (2011) 
Ponente et al. (2011) 

Singal et al. (2010) 
Padovani f20111 
Padovani QOll) 
Singal et al. (2010} 



^Calculated by integrating the number counts predicted by the iPadovani et al.l (1201 if ) 
observed SFG broken power law luminosity function at 1.4 GHz and assuming a = -0.8. 



'^Converted from 1 GHz using Eq. 19 and 20 of IPonente et al.l (120111 ). 



■^Converted from 1 GHz using Eq. 19 and 21 of iPonente et al.l (120111 ). 

"^Estimated based on the expectation that LRPEs have surface density and flux densities 
similar to radio quiet AGN. 



Note. — SF refers to star formation, SFG refers to star forming galaxy, IRC refers to the 
infrared-radio correlation, RSNe refers to radio supernovae, and LRPE refers to low radio 
power ellipticals. 
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Fig. 1. — Euclidean normalized 1.4 GHz differential number counts for bare AGN an d AGN 
with star f ormation. The blue lines show bare AGN assuming the iLa Franca et al.l (120051 ) 
HXLF and iPanessa et al.l (120071 ) Rx , as summarized in the fo urth line of Table [H Th e black 
lines sho w AGN with host SFR as de scribed by Equation^ the lLa Franca et al.l (120051 ) HXLF 
and the iTerashima fc WilsonI (120031 ) Rx , as summarized in the tenth row of T able [H Also 



shown are the observed source counts fr om a vari ety of surveys: VLA-C QSMOS (IBondi et al. 



20081). CDF-S (iKellermann et al.ll2008h . SSA 13 (iFomalont et al.ll2006h. 



( Sevmour et al.ll2C)04r) . Subur u/ XMM-Newton Dee p Field (ISimpson et al. 



3 hr XMM-Newton 



2006), VLA-VVDS 



( IBondi et al.ll2003[ ). Phoenix jHopkins et al.ll2003l ). and FIRST ( Iwhite et al.lll997l ). The blue 
and red points show the estimated AG N radio counts from t he 13 hr XMM-Newton (blue; 
Seymour et al. 2008) and CDF-S (red; Padovani et al. 20 111). The cyan point s show the 



estimated radio quiet AGN 1.4 GHz counts in the CDF-S (jPadovani et al.ll201l[ ). 
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Fig. 2. — Euclidean normalized differential number counts for the maximal AGN and host 
star formation contribution t o the 1.4 GHz CRB (last line of Tabled]). The radio loud AGN 
are accounted for using the ISmolcic et al.l (120091) PDE lu minosity function and the radio 
qu iet AGN are accounted for using the lUeda et al. I J2003h HX LF, t he evolving fcr model 
of braper fc Ballantvnel f |201oh . and the iTerashima fc WilsonI J^OOsI) Rx- The SFR of the 



radio quiet AGN hosts is described by Equation O This model gives the maximal total AGN 
and host brightness temperature of 0.042 K. The plotted data is the same as in Figure [TJ 



